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ABSTRACT 

We present a method to test the isotropy of the magnitude-redshift relation of Type la Supernovae 
(SNe la) and single out the most discrepant direction (in terms of the signal-to-noise ratio) with respect 
to the all-sky data. Our technique accounts for possible directional variations of the corrections for SNe 
la and yields all-sky maps of the best-fit cosmological parameters with arbitrary angular resolution. To 
show its potential, we apply our method to the recent Union2.1 compilation, building maps with three 
different angular resolutions. We use a Monte Carlo method to estimate the statistical significance 
with which we could reject the null hypothesis that the magnitude-redshift relation is isotropic based 
on the properties of the observed most discrepant directions. We find that, based on pure signal-to- 
noise arguments, the null hypothesis cannot be rejected at any meaningful confidence level. However, 
if we also consider that the strongest deviations in the Union2.1 sample closely align with the dipole 
temperature anisotropy of the cosmic microwave background, we find that the null hypothesis should 
be rejected at the 95 — 99 per cent confidence level, slightly depending on the angular resolution of 
the study. If this result is not due to a statistical fluke, it might either indicate that the SN data 
have not been cleaned from all possible systematics or even point towards new physics. We finally 
discuss future perspectives in the field for achieving larger and more uniform data sets that will vastly 
improve the quality of the results and optimally exploit our method. 

Keywords: cosmology:dark energy, supernovae:general, methodsidata analysis. 


1. INTRODUGTION 


In 1998, the luminosity-redshift relation (Hubble dia¬ 
gram) of a few tens of Type la supernovae (SNe) pro¬ 
vided the evi dence base for the accelerated expansion o f 


1998 Perlmutter et al. [1999 ) 


the universe (Riess et al. 

Since then, major efiorts have been made to increase the 
sample size, extend it to higher redshift, and refine the 
observational and data-reduction techniques. Current 
datasets already include several hundreds of objects but 
the quest for dark energy drives copious activity in this 
field. 

The control of systematic errors is the key to making 


the study o f SNe la a prime cosmological tool. [Suzuki 


et al. (2012) state that systematic uncertainties already 


dominate over the statistical ones in the determination 
of the cosmological parameters. Given that systematics 
will become even more important in the future, a care¬ 
ful scrutiny of all the possible sources of methodological 
bias is crucial. In this paper, we focus on the spatial 
isotropy of the Hubble diagram traced by type la SNe. 
The standard cosmological model is rooted in the as- 
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sumption that the Universe is homogeneous and isotropic 
on large scales. Hence, SNe la are expected to (statisti¬ 
cally) obey the same dimming relation in all directions. 
There are, however, several phenomena that could intro¬ 
duce anisotropies with different characteristic scales and 
amplitudes in the observed expansion rate. To name a 
few: dust absorption (both in the Milky Way and in the 
galaxies hosting the SNe), redshift-space distortions due 
to large-scale motions, weak gravitational lensing, the 
presence of large-scale structures and contamination of 
the SNe la samples. Detecting these effects and correct¬ 
ing for them would ultimately lead to tighter and less 
biased constraints on the cosmological parameters. 

At the same time, it is healthy to scruti nise the valid¬ 


ity of the standard model of cosmology ([Kroupa 2012 


__ __ dK 

Kroupa et al.[[2Q12 Kroupa[[2Q15 Koyama|[2015 ) and its 


fundamental assumptions, namely those of the cosmo¬ 
logical principle. Ruling out cosmic isotropy with high 
statistical confidence would lead to a major paradigm 
shift especially if such a conclusion is confirmed by mul¬ 
tiple datasets affected by different systematics. In this 
respect, the analysis of temperature anisotropies in the 
Cosmic Microwave Background (CMB) has dominated 
the scene in the last decade. The WMAP satellite de¬ 
tected a few large-scale “anomalies” that somewhat de¬ 
viate from the expectations of the standard model that 
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best fits the data on smaller scales (|Teg:mark et al. 2003 


Eriksen et al. 2004 Hansen et al.|2009| Copi et al.|2010a| 

in brief, the quadmpole and octopole terms are surpris¬ 

ingly planar and there is a significant alignment between 
them. Moreover, their normals he close to the axis of 
the CMB dipole. This discovery generated a long lasting 
debate in the literature concerning whether or not these 
features are genuine signs of new physics. 

Alternatively they could be due to the influence of 
data processing, to the imperfect removal of foregr ound 
contaminants a nd secondary astrophysical effects (|Ras- 
sat et al. 12014), as well as to a statistical fluke (IBen- 


nett et al. 


. The Planck satel lite recently confirmed 


the existenc e of these alignments (Planck Collaboration 
et al.|2014b|) suggesting that they are not artifacts of the 
data-reduction pipelines. A satisfactory explanation for 
the origin of these asymmetries is still not available. 

The isotropy of the H ubble diagram for SNe l a has 
been repeate dly tested. Kolatt & L ahav ( '2001|) used 

( 1998| ) and Perlmntter et al. 


79 SNe from Riess et al. 


(1999) to perform localized fits within an opening angle 
of 60"^ around random directions. After expanding the 
best-fitting cosmological parameters in low multipoles, 
they found that no dipole anisotropy was statistically sig¬ 
nificant. Subsequent studies mainly adopted two meth¬ 
ods: either they compared Hubble diagrams for pairs of 
hemispheres and looked for the most discrepa nt hemi¬ 
spheric cut (Hemi spherical Comparison, e.g. Schwarz 
& Weinhorst||2007|) or fit a dipo le angular distribution 
(Dipo le Modulation Fitting, e.g. Cooke & Lynden-Bell 
2010). Other au thors looked for angula r correlations in 
SN rnagnitudes (Blomqvist et al. 2008) or analyzed the 
magnitnde-redshnt re lation in the context of anisotropic 
cosmological mo dels (Koivisto & Mota 2008 Campan- 
elli et al. 2011). Low-redshift samples were used to 
estim ate the direction and amplitude of the local bulk 
flow dBonvin et al. 1 120061 iSchwarz & Weinhorstl 12007 


loo 


Colin et al.| 

pnnr 

'Turnbull et al.| 2U12| |Rathaus et al. 

20131 IFeindt et al. 

20131 Kalus et al. 2013 Appleby & 


several authors analysed h^ 


ler-redshift data to look 


for large-scale anisot ropies dSchwarz fc WeinhQr^ 2007 


Gupta & Saini |20101 

Cooke ^ Lynden-Bell||20101 fAnto- 

niou & Perivolaropou 

os|20101 Mariano & Perivolaropou- 


2011[ |Zhao et al.||2013[ [Heneka ~ 


Wang|20T^ Yang et al.|2014f Chang fc Lin|2015 Jimenez 


.1120141 IWang & 


et al. |2015|), which is also the aim of our work. Statisti¬ 

cally signmcant deviations have b een detected at low red- 
shift (Schwarz & Weinhorst|2007), while no high-redshift 
study could rule out isotropy at more than 2 Gaussian 
standard deviations, cr. 

In this paper, we present a simple but powerful method 
to test the isotropy of the luminosity-redshift relation for 
SNe la. Contrary to most previous studies, our analysis 
neither searches for hemispheric asymmetries and dipolar 
patterns nor does it use any other template anisotropic 
configuration. For each direction on the celestial sphere 
r G we derive a set of cosmological parameters by 
considering only the SNe that lie within an angle 0 from 
r. We then build maps of these “local cosmological pa¬ 
rameters” with different values of 0 and identify the di¬ 
rections associated with the most significant anisotropies 
taking into account that the number of datapoints used 
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Figure 1. Mollweide projection map of the Union2.1 SNe la in 
the Galactic coordinate system. Each circle corresponds to the 
position of a SN la on the sky and is colour-coded based on the SN 
la redshift. The black solid curve indicates the celestial equator. 

in the fit fluctuates from one direction to another. For 
completeness, we consider that the correction for the dis¬ 
tance modulus of SNe la might also depend on r due, for 
instance, to dust extinction. Therefore, our strategy is 
able to detect anisotropies generated both by physical 
effects and by systematics. Even though our method is 
best suited for the large SN samples with nearly uniform 
sky distribution that will become available in the next 
decade, we provide an example of its poten tial by apply¬ 
ing it to the Union2.1 SN la compilation (Suzuki et al. 
2012) from the Supernova Cosmology Project (SCF). We 
limit our study to redshifts ^ > 0.2 in order to minimize 
the influence of local inhomogeneities and bulk flows. 

The rest of this paper is organised as follows. Section 
2 describes the main properties of the Union2.1 sample. 
Our method of analysis is introduced in Section 3. Re¬ 
sults are presented and critically discussed in Section 4. 
Finally, we conclude in Section 5. 

2. DATA 

The Union2.1 compilation (Suzuki et al.||2012) collects 
data for 580 SNe la in the redshift range of 0.015 < 
2 : < 1.414. It combines entries from 19 datasets uni¬ 
formly analysed after adopting strict light curve quality 
cuts and the SALT2 lightcurve-fitter (|Guy et al.||2007|. 
The Union2.1 catalog has been built for dark-ene rgy sci¬ 
ence and upda tes the previ ously released Union ([K owal¬ 
ski et al.|2008|) and Union2 (jArnanullah et al.|2010l 
pilations. In particular, it contains 14 new SlNe 


com- 
iscov- 

ered in the HST Cluster Supernova Survey (a survey run 
by the SCP) that pass the Union2 selection cuts. Ten of 
these SNe are at z > 1 which makes the Union2.1 sam¬ 
ple ideal for studying isotropy out to the largest possible 
distances. 

The Union2.1 catalog provides five entries for each 
SN, specifically: name, redshift (CMB centric), distance 
modulus, error in the estimate of the distance modulus, 
and the probability that the SN was hosted by a low- 
mass galaxy. We obtained the coordinates for all the 
SNe la in the compilation either from the NASA/IPAC 
Extragalactic Database (NED^or directly from the Su- 

^ The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech¬ 
nology, under contract with the National Aeronautics and Space 
Administration. 









































































































































































Probing isotropy with SNe Ia 


3 



Figure 2. Redshift distribution of Union2.1 SNe. 


per N ova Legacy Survey (SNLS) data release (Astier et al. 
2006|). The sky distribution of Union2.1 SNe is plotted in 
b'igure The angular position of each SN is marked by 
a symbm which has been colour-coded based on redshift. 
Several features are immediately apparent in the image. 
First, there are only a few SNe close to the galactic plane. 
Second, an arc-like region in the southern hemisphere 
is much more densely populated than the rest. This is 
the footprint of the Sloan Digital Sky Survey-II (SDSS- 
II) SN search and corresponds to the southern equatio- 
rial stripe (Stripe 82, with coordinates — 50 < RA < 59 
and —1.25 < DEC < 1.25) which has been imaged re¬ 
peatedly with broad wavelength coverage and also been 
subje ct to extensive spectroscopic studies (Kessler et al. 
2009). Finally, high-redshift SNe are very sparsely dis¬ 
tributed and rare which is also evident from the redshift 
distribution of the Union2.1 SNe shown in Figure 


3. METHOD 

As mentioned in the Introduction, previous studies on 
the isotropy of the luminosity-redshift relation of SNe la 
have mainly searched for dipolar anisotropies or hemi¬ 
spheric asymmetries. In this section, we introduce a more 
general method that does not assume any particular form 
of the anisotropy. 


3.1. Cone Analysis 

Let us consider a particular direction on the sky, f, 
with Galactic longitude I and latitude b. In order to 
single out a finite region surrounding f, we consider a 
cone with apex angle 20 subtending a solid angle of 
^cone = 27 r(l — COS 6>) sr on the celestial sphere. The 
apex of the cone is located at the centre of the Galac¬ 
tic coordinate system and its axis of symmetry points 
towards f. After isolating the SNe la contained within 
the cone, we build their magnitude-redshift relation and 
derive “local cosmological parameter s” b y fitting a theo¬ 
retical relationship to it (see Section 3.2 for details). We 
then vary the cone direction r making sure that we cover 
the whole sky. For convenience, we move r along t he 
pixel centers of a HEALPi:x[^grid ( |G6rski et al.||2QQ5 ). 


^ Hierarchical Equal Area isoLatitude Pixelization, 
http://healpix.sourceforge.net 


For the application of the method to the Union2.1 data, 
we repeat the analysis using three different opening an¬ 
gles: ^ = f (hemispheres), | and We use a HEALPix 
grid with 192 pixels so that the solid angle subtended by 
each pixel is much smaller than that subtended by the 
cones. 


3.2. Formulation 
3.2.1. Global fit 

SNe la are not perfect standard candles, their peak 
brightness correlates with their color, the light-curve 
width and the mass of the host galaxy. In the Union2.1 
sample, individual lightcurves are analyzed with the 
SALT2 fitter which provides estimates for three param¬ 
eters: the peak magnitude, niobs^ in the rest-frame B 
band, the deviation, xi, from the average light-curve 
shape and the deviation, c, from the mean B — V color. 
The color and light-curve-shape corrected distance mod¬ 
ulus is then written in terms of four unknown parameters 
and Mb) so that 

fJjB ex ' X\ P ' C S ' Phost Af^ , (l) 

where Mb is the absolute B-band magnitude at max¬ 
imum of a SN la and Phost denotes the probability 
that the SN la is hosted by a galaxy with stellar mass 
< IO^^Mq. This probability is estimated differently 
for untargeted and targeted surveys. 

In the context of the theory of general relativ¬ 
ity, homogeneous and isotropic universes are described 
by Eriedmann-Lemaitre-Robertson-Walker models. Eor 
simplicity we only consider flat models in which the den¬ 
sity parameters for the matter and the cosmological con¬ 
stant satisfy the relation + ^a = 1- Eollowing stan¬ 
dard practice, we write the magnitude-redshift relation 
of SNe la in terms of the distance modulus 


li{z) = 5 logio dUz, Oa) + 5 logio ’ (2) 

where 


dr = 


/■ 


(1 -h z) dq 


(1 -j- g) (1 + ^ttiQ) — ^(2 + q)^A 


nl/2 


(3) 


is the dimensionless “Hubble-constant-free” luminosity 
distance and Dh = c/Hq is the Hubble radius defined in 
terms of the speed of light and the present-day value of 
the Hubble constant. 

Classically, the SN data are fitted with a cosmological 
model assuming Gaussian err ors and following a maxi¬ 
mum likelihood approach (e.g.|Astier et al.|2006|). Eor N 
Type la SNe, this corresponds to minimising the target 
function 

^2 = C-1V (4) 


where V is a Wdimensional vector with elements Vi = 
IXBh{(x^ P^S^Mb) — and C is the covari¬ 

ance matrix of the errors in the observed distance mod¬ 
uli. Eor the Union compilations, this matrix is pub¬ 
licly available. Its off-diagonal elements include sev¬ 
eral contributions due to the light-curve fits, galactic 
extinction, gravitational leasing, peculiar velocities and 
sample-dependent systematics. The nuisance parameters 
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a, /3,S and Mb are fitted simultaneously with the cosmo¬ 
logical parameters. Actually, the best-fit values for Mb 
and Hq are completely degenerate as only the combina¬ 
tion M = Mb + 5\ogiQ{DH/Mpc) appears in eq. 0. 
Usi ng the whole data s et gives the following best-fit 
ues ( Suzuki et al.|2Ql^ a = 0.121, P = 2.47, 5 = —0.032, 
and Mb = —19.321 (for Mq = 70 km s“^Mpc“^). 


3.2.2. Local fits 

The parameters a, (3 and S describe correlations be¬ 
tween different SN observables a nd might vary for the 
different surveys of a compilation. Karpenka et al. (2015) 
found inconsistencies between the values of these correc- 
tion parameters in the Union2 catalog. In order to ac¬ 
count for possible direction-dependent systematics, when 
we consider localized sub-samples of the Union2.1 data, 
we should in principle allow them to vary freely. How¬ 
ever, this would require knowledge of the covariance be¬ 
tween rriobs^ and c for individual SNe. Regrettably 
this information is not provided in the Union2.1 catalog. 
We therefore adopt a simplified approach by assuming a 
constant correction, jHcor^ for fhe distant modulus of all 
the SNe lying within a cone. In other words, we keep 
the quantities a, 6 and A4 fixed at their global best-fit 
value (hereafter denoted with a hat) but we write 


Vi = /?, S, Mb) - 5 logn, diizi; Qa) ~ Aq (5) 

with Aq = 5 log]^Q(ciJ|7^) + 25 — Hcor- Note that Aq ac- 
counts for both an “anisotropic Hubble constant” and 
for the mean effect of variations in and Mb due 

to systematic errors. We are left with a two-dimensional 
problem. For each pixel on the sky, we then determine 
the best-fitting values of the cosmological parameter Ua 
and of the correction parameter Aq by minimizing the 
target function (covariances are extracted from C after 
identifying the SNe in the cone). However, the model pa¬ 
rameters anticorrelate: directions associated with large 
values of Oa provide low values of Aq (and viceversa). In 
order to minimise this effect, we use the luminosity dis¬ 
tance, evaluated at the mean redshift of the sample 
{z = 0.36) as a pivot point and define 

Vi = - Slogjo ~ ^ 

where A = Aq + 5 log;L0 ^e(^; ^a) ^a are our free 
parameters. 


4. RESULTS 

4.1. All-sky fit 


To test the consistency of our approach with previous 
studies, we first perform an all-sky fit. Results are shown 
in Table[^for the entire Union2.1 sample and for two sub¬ 
sets including the SNe la with redshift smaller and larger 
than 2 : = 0.2 (in this paper, uncertainties on the value of 
single parameters always correspond to Ay^ = !)• Our 
results are in excell ent agreement with the analysis in 
Suzuki et al. (2012) who found Ua = 0-705tQ Q 43 
their Table 7). Also note that setting i^cor = 0, Mq = 70 
kms“^Mpc“^ and Mb = —19.321 corresponds to Aq = 
43.159 mag which gives A = 41.419 mag for Ua = 0.705 
and z = 0.36. 


Table 1 

All-sky fitting results 




A (mag) 


All SNe 

^ > 0.2 

2 ; < 0.2 

''^-0.04 

n fisT0.06 

41.428+"o^"3« 

41.380t°;MO 

0.94 

0.94 

0.93 


Note. — Best-fit parameters and the corresponding reduced chi- 
square, tor the entire Union2.1 sample and for two redshift 

subsets. The quoted uncertainties correspond to Ax^ = 1- 


4.2. Cone analysis 
4.2.1. Qa maps 

Sky maps of the best-fit values for Ua (left) and A 
(right) are shown in Figure for three different cone 
opening angles (from top to bottom: d = |, | and | ra¬ 
dians). These have been obtained using all the Union2.1 
SNe with redshift z > 0.2. White pixels indicate the 
directions (mostly located around the Galactic equator) 
in which the corresponding cone contains less than 25 
SNe la. These directions are excluded from all statisti¬ 
cal analyses because they are associated with extremely 
large errors in the fitted parameters. Of course their 
number increases with decreasing the opening angle of 
the sampling cone. Similarly, the size of fiuctuations in 
the best-fit values for Ua and A increases with reducing 
0 . 


4.2.2. Most discrepant directions 

Although Figure gives a first visual impression of the 
local best-fit parameters, it does not take into account 
the non-uniform sky coverage of the Union2.1 data set. 
For a given opening angle, different directions on the ce¬ 
lestial sphere are generally associated with very different 
numbers of SNe la. This strongly influences the uncer¬ 
tainty of the best-fit values. 

In order to single out the most discrepant directions 
in a statistically meaningful way, we assume the null hy¬ 
pothesis that the Universe follows the cosmological prin¬ 
ciple and there are no angle-dependent systematic effects 
plaguing the Union2.1 sample. For each pixel we then 
evaluate the target function fixing the free parameters 
at the values Cl a = 0.70 and A = 41.428 that provide the 
best-fit solution for the complete Union2.1 sample. How¬ 
ever, only the SNe within the sampling cone are used to 
calculate the value that we denote by Finally, we 
estimate the probability P that random noise could gen¬ 
erate a value exceeding Assuming Gaussian errors, 
this probability coincides with the cumulative chi-square 
distribution function evaluated at 


P = 


1 

2fo2r(zy/2) 


POO 


(7) 


where v is the number of degrees of freedom - i.e. the 
number of SNe la used in the fitting procedure minus two 
(the number of free parameters). We adopt the P value 
as a measure of how well the all-SNe best-fit parameters 
also describe the SN data in a specific direction on the 
sky. Gonsequently we identify the most discrepant di¬ 
rection (i.e. the direction showing the most statistically 
significant deviation from isotropy) with the pixel show¬ 
ing the smallest P value. It is worth stressing that this is 
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Figure 3. Best-fit values for Oa (left) and the correction parameter A (right) in different directions on the sky. Each pixel shows results 
that have been determined considering the magnitude-redshift relation of all Union-2.1 supernovae with redshift 2 : > 0.2 that he within 
an angle 6 from the pixel center. The cone opening angle 6 assumes the values ^ (top), ^ (middle) and ^ (bottom). The white regions 
indicate the directions for which the sampling cone contains less than 25 SNe la. 


not necessarily the direction in which the Universe (or the 
SN data) might present the strongest intrinsic anisotropy 
but only the direction in which, given the current data, 
we can measure the most meaningful deviation in terms 
of the signal-to-noise ratio. 

Maps of the P value are plotted in Fig. [^for the three 
different cone opening angles. The most (second-most) 
discrepant directions are highlighted with a star (circle) 
in each panel. Further information is provided in Table 
1^ which gives the P value, the coordinates and the num- 
b^er of SNe la associated with the most and the second- 
most discrepant directions together with the local best-fit 
parameters. The motivation for showing two directions 
per map is as follows: i) the difference between their P- 
values is small (see Table |^, ii) the covariance matrix 


provided by the SCP is likely to be a noisy estimate, and 
iii) neglecting off-diagonal covariances switches the order 
between these directions for ^ 

Intriguingly, the most discrepant directions obtained 
with the three cone opening angles lie close to each other. 
Also the best-fit parameters are quite similar (let us not 
forget, however, that the maps with different 0 are not 
independent as they use the same SNe and that there is 
significant overlap between the most discrepant cones). 
In Figurej^we compare the formaQla confidence regions 
(Ay^ < 2.30) obtained from the all-SNe fit against those 
derived from the local fits along the most discrepant di- 

I.e. obtained assuming independent Gaussian errors . The lim¬ 
itation of this approach is discussed in detail in Section |433| 
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Table 2 

Most discrepant directions 


0 

{i,b) 


A 


P 

N 

★ 7r/2 

(33.7,-19.5) 


41 

^-^•^^^-0.072 

1.09 

0.192 

128 

• 7r/2 

(0.0,-30.0) 

0 . 611 Ho 

41 441+0-066 

1.08 

0.197 

161 

★ 7r/3 

(112.5,-9.6) 

0.60+^;^^ 

41 43 Q+o-J^ 6 u 

1.20 

0.086 

74 

• 7r/3 

(56.2,-41.8) 


41 424+0-081 
ni.nzn Q Q72 

1.15 

0.101 

118 

★ 7r/6 

(67.5,-66.4) 


41 419+0-084 

1.18 

0.081 

100 

• 7r/6 

(101.2,-41.8) 

0.55l“Ji 

41 408“^®'^^'^ 

1.20 

0.085 

73 


Note. — Galactic coordinates (/, b) and P values characterizing 
the most (stars) and the second-most (circles) discrepant directions 
for different cone opening angles, 6. Also reported are the number 
of SNe la in the cones, A/”, the best-fit values for and A (in 
mag) and the ratio x^/^- 



Table 3 

Results of NGH and SGH fitting 



11a 

A 

xV^ 

N 


NGH 

SGH 

0 70+O-U7 
O'- ' 0>_o.09 

0'-0>O 0.08 

41 44Q+U.U88 

ni.nnc)_Q Qg4 

41 441+0.060 

ni.nni_Q Qgg 

0.82 

1.01 

123 

227 

0.30 

0.23 


Note. — Best-fit values obtained using SNe in the NGH and 
SGH, separately. Also reported are the corresponding reduced chi- 
square, x^/^: fhe number of SNe la with 2; > 0.2 considered for 
the fit, A, and their average distance-modulus uncertainty, 


rections. In all cases, the tension between the local and 
the global fits is marginal and the formal la regions al¬ 
ways overlap. 

Visual inspection of Figure shows a striking con¬ 
trast between the P values measured in the Northern and 
the Southern Galactic Hemispheres (hereafter NGH and 
SGH, respectively), although there is no tension between 
the luminosity-distance relation in the two hemispheres 
(see Table . The discrepancy in the P values is mainly 
due to th^act that the Union2.1 uncertainties in the 
distance-moduli are on average 30 per cent larger in the 
NGH. Consequently, the reduced chi-square /i' tends 
to be smaller in the NGH although there are many more 
SNe in the SGH (227 vs 123) to drive the fit results for 
SNe with 2 ; > 0.2 closer to the SGH results. 

4.2.3. Monte Carlo analysis 

Taken at face value, the probabilities P associated 
with most discrepant directions (see Table are moder¬ 
ately significant. However, assuming Gaussian errors is 
a strong limiting factor. Also, the size of the errorbars in 
the distance modulus (and the off-diagonal covariances) 
provided in the Union2.1 catalog might be inaccurate 
and, as a consequence, inference based on the statistic 
might be biased. For these reasons, we re-evaluate the 
statistical significance of the anisotropies using a more 
robust Monte Carlo method. 

In order to assess the impact of random errors and 
account for the non-uniform angular distribution of the 
Union2.1 sample, we build 1000 mock catalogs by ran¬ 
domly shuffling the distance moduli of the Union2.1 SNe. 
In practice, we assign the distance modulus, its uncer¬ 
tainty and the redshift of a SN la to the angular position 
of another (random) SN la. Each mock catalog thus 
contains exactly the same number of SNe as the origi¬ 
nal Union2.1 sample and has exactly the same SN sky 




0 1 

p 

Figure 4. Maps of the P value (the estimated likelihood of getting 
larger deviations than in the data due to random fiuctuations under 
the null hypothesis that the SN la Hubble diagram is isotropic) for 
different cone opening angles (from top to bottom, ^ = f 5 f and 
^). The direction with the smallest value of P in each map is 
marked with a star and the second-most discrepant direction is 
highlighted with a circle. The white regions denote the pixels for 
which the cone contains less than 25 SNe la and are excluded from 
the statistical analysis. 


distribution. Moreover, all possible anisotropies should 
be erased by the shuffling procedure while the statistical 
properties of the distance moduli and their uncertainties 
are unchanged with respect to the observational data. 
Therefore our mock catalogs form an ensemble of real¬ 
izations mimicking an isotropic Universe but having the 
same statistical properties as the actual Union2.1 data. 

We treat the mock catalogs as the real data and iden¬ 
tify the two most discrepant directions in each of them 
using the P-value-based method for the three different 
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cone opening angles. We then compute the fraction, 
/o, of the realizations in which the most (or the sec¬ 
ond most) discrepant direction is associated with a P 
value which is smaller than the observed ones reported 
in Table For the most (second-most) discrepant di¬ 
rections we find that /o = 0.569, 0.623 and 0.674 (0.473, 
0.550, 0.490) for 0 = | and respectively. Purely 

based on this signal-to-noise criterion, we conclude that 
no statistically significant anisotropy can be detected in 
the Union2.1 sample. 


4.2.4. Alignment with the CMB dipole 

Although the Monte Carlo test shows that random 
chance in an isotropic universe can easily produce most 
discrepant directions with lower P values than we found 
analyzing the actual data, the observed anisotropies 
present a characteristic feature which is worth being dis¬ 
cussed. 

The temperature distribution in the CMB presents 
a strong dipole anisotropy which is usually interpreted 
as due to our motion with respect to the CMB rest 
frame towards the direction with Galactic coordinates 
(263°9 9 T 0.14,48?26 T 0.03) ( [Planck Collaboration et aT 


2014a). Figureshows that the most discrepant direc- 


tions we obtain^ from the Union2.1 sample closely align 
with the axis of the CMB dipole (CDP) in the SGH op¬ 
posite to our motion with respect to the CMB rest frame 
(hereafter CDP-South). Assuming that the redshifts of 
the SNe la in the Union2.1 compilation have been cor¬ 
rectly transformed to the CMB rest frame, there is no 
obvious reason for explaining the origin of this align¬ 
ment. As already mentioned in the Introduction, the 
CMB quadrupole (CQP) an d octopole (COP) are also 


closely aligned with the CDP ([Planck Collaboration et al. 
2Q14b[ l^chwarz et al.[[2QQ4| |Copi et al.||2Ul(Jb[ |2(J13[ see 
Figure!^. It is yet unclear whether these alignments are 
a statistical fiuke or a signature of new physics. Any¬ 
way, our study shows that the magnitude-redshift re¬ 
lation of SNe la with 2 ; > 0.2 tends to be different in 
the same direction (albeit the difference is detected with 
low signal-to-noise ratio). Other authors have reported 
similar results usi ng the Union compilations ([Cooke fc 


Lynden-Bell|20I01|Antoniou & Perivolaropoulos|2UIUt|Li 

et al.[[20I3[). 

4'he debate on the physical relevance of the CMB 
anomalies opened up a discussion in the literature about 
the legitimacy and validity of “a posteriori” analyses in 
which tailored statistical tests are designed and hand 
picked after noticing the peculiarities in the data. A 
widespread point of view states that in a large dataset it 
is always possible to isolate some “strange” features (e.g. 
Bennett et al.[[20II[). To minimize the pitfalls of a pos¬ 
teriori reasoning, we focus on the well established CDP 
and do not consider the CQP and COP any further. 

We thus proceed to quantify the probability that the 
most discrepant directions (defined in terms of the P 
value as above) form a given angle with the CDP under 
the null hypothesis of an isotropic magnitude-redshift re¬ 
lation. In order to account for the non-uniform sky distri¬ 
bution of the Union2.I sample (especially for the SDSS-II 
stripe which is close to the CDP-South) we use the Monte 
Carlo realizations introduced in Section [4.2.3[ Figure 
shows the resulting probability distribution for the cosine 
of the angle between the most discrepant direction and 


the axis of the CDP-South. Our measurements from the 
Union2.I data are indicated by vertical dashed lines. The 
second column in Table [^ reports the fraction of Monte 
Carlo realisations, /i, showing a better alignment than 
our measurement. Our results suggest that it is rather 
unlikely to get an alignment as strong as the observed 
one under the null hypothesis of an isotropic magnitude- 
redshift relation. In fact, considering the most discrepant 
direction for 6> = only 8.5 per cent of the Monte 
Carlo realisations show a smaller separation angle than 
observed and this reduces to 4.5 percent for ^ Note 
that for 0 = ^ and the second-most discrepant direc¬ 
tions are even better aligned with CDP-South. In these 
cases /i ^ O.OI. 

The test above is blind to the statistical significance 
of the most discrepant directions. In order to account 
for this, we compute the fraction of Monte Carlo reali¬ 
sations, / 2 , for which the most discrepant directions are 
at least as significant as the measured ones (in terms 
of the P value) and are also better aligned with the 
CDP-South. The third column in Table [4] shows that for 
the most discrepant directions this probability is smaller 
than 4.5 percent for all the opening angles which means 
the null hypothesis of an isotropic magnitude-redshift re¬ 
lation should be rejected at the 95 per cent confidence 
level. The value of /2 reduces to a fraction of a percent 
when considering the second-most discrepant direction 
for 6> = I and |. 

The measured anisotropy could be due to a statistical 
fiuke, to systematics in the SNe data (or error bars), to 
the presence of localized large scale structures, or even 
a sign of the failure of the cosmological principle. To 
further investigate its properties, we repeat the analysis 
along the most discrepant directions after slicing the SNe 
data in five redshift bins (0.2 < 2 ; < 0.3, 0.3 < 2 ; < 0.4, 
0.4 < 2 : < 0.6, 0.6 < 2 : < 0.9 and 2 : > 0.9). Regret¬ 
tably, due to the low number of SNe in each bin, the for¬ 
mal Icr errors span most, if not all, the parameter space 
0 < Da < I. Therefore no meaningful statements can be 
made regarding the variations of the best-fit cosmologi¬ 
cal parameters along the most-discrepant directions. In 
terms of signal-to-noise ratio, however, the redshift range 
0.4 < ^ < 0.6 clearly emerges as the most discrepant one 
for all the cone opening angles (P = 0.033,0.021 and 
0.016 for 0 = 7r/2,7r/3 and 7r/6, respectively). 

Given the current sparsity of the data, no firm con¬ 
clusion can be drawn except from the fact that there 
seems to be a moderately statistically significant (2-3cr) 
anisotropy in the magnitude-redshift relation of SNe la 
close to the direction opposite to our motion with respect 
to the CMB rest frame. It is worth remembering that, 
in the Union2.I compilation, most of the SNe la sur¬ 
rounding the CDP-South come from the SDSS-II stripe. 
Further investigations are thus needed to clarify the rela¬ 
tion between the CMB dipole axis, our motion, and the 
way SNe data around this direction are treated. 

5. CONCLUSIONS AND FUTURE PERSPECTIVES 

We presented a simple but powerful method for investi¬ 
gating the isotropy of cosmic acceleration traced by Type 
la SNe with different angular resolution, 0. The key idea 
is to consider all the SNe contained within a cone with 
vertex located at the origin of the Galactic coordinate 
system and with apex angle 20. “Local cosmological pa- 
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Figure 5. Formal Icr (Ax^ < 2.30) confidence regions from the all-SNe fit (gray) and the localized fits along the most discrepant directions 
(orange). From left to right the cone opening angle assumes the values ^ ^ and 



Figure 6. The most (star) and the second-most (circle) dis¬ 
crepant directions in the magnitude-redshift relation of SNe la 
(for three cone opening angles ^ , f and ^) obtained in this 

work are compared with the directions of th e CMB dipole fCDPb 
quadrupole fC QP) and octopole (COP) from|Planck Collaboration! 
|et al.|p014b| >. The black solid curve denotes the celestial equator. 


Table 4 

Alignment with the CMB dipole 



0 

OL 

h 

/2 

★ 

~T[2 

49°4 

0.085 

0.027 

• 

■k/2 

64°3 

0.152 

0.052 

★ 

7r/3 

45°5 

0.045 

0.021 

• 

7r/3 

20°5 

0.008 

0.002 

★ 

7r/6 

20"1 

0.064 

0.045 

• 

7r/6 

13°7 

0.010 

0.006 


Note. — Angular separation, o, between the direction of the 
CMB dipole in the SGH (CDP-South) and the most (stars) and the 
second-most (circles) discrepant directions for the maps based on 
the Union2.1 data with different cone opening angles, 6. The prob¬ 
ability of measuring a value smaller than a in random realisations 
of a isotropic magnitude-redshift relation is indicated with fi while 
/2 also accounts for the condition that the most (second-most) dis¬ 
crepant direction is associated with a smaller P value than for the 
Union2.1 measurement. Both probabilities have been estimated 
with a Monte Carlo method (see the main text for the details). 


rameters” are derived by fitting the magnitude-redshift 
relation of the SNe in the cone with a theoretical relation. 
The cone direction is then changed so that to cover the 
entire sky. Our cone-analysis method takes into account 
the mean variation of the SNe la correction parameters 
over different directions, and yields all-sky maps of the 
best-fit cosmological parameters. 


Although a large data set with a uniform sky distribu¬ 
tion is required for a thorough investigation of isotropy, 
we provided an example of the potential of our method 
by applying it to the SNe la with redshift ^ > 0.2 in the 
Union2.1 compilation. Assuming a flat Universe in the 
context of the standard cosmological model, we fitted 
the magnitude-redshift relation by varying the density 
parameter of the cosmological constant, Ua, and a pa¬ 
rameter , A, including the effect of both the Hubble con¬ 
stant and the mean SNe la correction parameters. We 
used a HEALPix grid to discretise the celestial sphere 
and obtained sky maps for Ua and A considering three 
different cone-opening angles ^ = f, f and 

We ranked the pixels in each map in terms of a P 
value derived from the distribution and which mea¬ 
sures how much the local fits differ from the cosmology 
determined using the entire Union2.1 sample (in a signal- 
to-noise sense). We thus found the most discrepant di¬ 
rections (two per cone opening angle). Finally, we used 
a Monte Carlo method to estimate the statistical sig¬ 
nificance at which we could reject the null hypothesis 
that the magnitude-redshift relation of SNe la is isotropic 
based on the properties of the most discrepant directions. 
We found that random fluctuations can easily produce 
deviations from isotropy with smaller P values than mea¬ 
sured in the Union2.1 data. Therefore, the null hypothe¬ 
sis cannot be rejected at any meaningful confidence level 
based on signal-to-noise arguments alone. However, if 
we also consider that the detected anisotropies in the 
Union2.1 sample align well with CMB dipole axis in the 
Southern Galactic Hemisphere, we find that the null hy¬ 
pothesis should be rejected at the 97.3, 97.9 and 95.5 per 
cent confidence level for opening angles ^ = f, f and 
respectively. 

We conclude that, although the deviation from 
isotropy that we found is not very significant per se in 
terms of signal-to-noise ratio, its vicinity to the axis of 
the CMB dipole (which enters the pipeline to determine 
the SN redshift in the CMB rest frame) with 2-3cr statis¬ 
tical significance requires further investigation both on 
the observational and on the theoretical sides. Note that 
other observations detected anisotropies in the same area 
of the sky. The statistical significance of the quadrupole- 
octopole alignment in the CMB is approximately 99 per 
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cos(a) 
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Figure 7. Top: Probability distribution of the cosine of the angle between the direction of the CMB dipole in the SGH, rcDP-South^ cind 
the m ost discrepant direction of the Hubble diagram of SNe la, Vdis: determined using the Monte Carlo realisations introduced in Section 
|4.2.3| Bottom: As above but with the additional condition that the most-discrepant direction is associated with a smaller P value than 
for the Union2.1 measurement. From left to right, the panels refer to the cone-opening angles ^ J and The cosine of the observed 

angular separation in the Union2.1 sample, a, (reported in Table is shown as a dashed line. 



Figure 8. Survey footprints for the Euclid deep field, LSST main 
survey and the (likely) SNe fields of DES in the Galactic coordinate 
system. The size of the fields for DES and Euclid have been artifi¬ 
cially magnified to ease readability. The black solid curve denotes 
the celestial equator. 


cent (Planck Collaboration et al. 2014b). On combina¬ 
tion of the likelihoods between the CMB and SN la, the 
null hypothesis of isotropy should be rejected at the 99.98 
per cent confidence level (approximately 3.5 Gaussian a). 
In this paper we followed a conservative approach by only 
considering the SN la data. 

This study should be repeated when larger data sets 
with more uniform sky coverage will be available. Sev¬ 


eral major current and future facilities have dedicated 
plans for studying the accelerated expansion of the uni¬ 
verse using SNe la. For instance, the Dark Energy Survey 
(DES) integrates a dedicated program that should detect 
around 4000 SNe la in the redshift range 0.05 < ^ < 1.2 
(jBernstein et al.||2012). Similarly, the Euclid mission in- 
cludes a SNe survey within two deep fields each covering 
around 20 deg^ and is expected to discover about 3000 
SNe la out to z ^ 1.2 (Laureijs et al. 2011). However, 
both these surveys will only provide SlN data in rela¬ 
tively small regions of the sky (see Eigurel^ and the most 
promising perspective for isotropy tests oi ]he Hubble di¬ 
agram comes from the Large Synoptic Survey Telescope 
(LSST). While its use for a SN-dedicated survey on a lim¬ 
ited area of sky will be able to deliver as many as 140, 000 
SNe la (in 10 years) with very precisely measured light 
curves, in its normal operating mode (due to its rapid 
cadence), LSST will discover around 250, 000 SNe la per 
year in the redshift range 0.45 < z < OT a nd across 
a large fraction of the sky (Ivezic et al. [2008). Einally, 
the Panoramic Surv ey Telescope ^ K apid Response Sys¬ 
tem (Pan-STARRS,|Brout et al.|2013|) which is observing 
the Northern part of the sky will complement the above 
mentioned surveys. In summary, exciting perspectives 
to test the isotropy of the magnitude-redshift relation of 
SNe la with unprecedented accuracy will open up within 
the next two decades. 
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